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ABSTRACT 

^  Six  degree-of-freedcm,  rigid  body  equations  of  motion  are  described  suitable  for 
modeling  the  dynamic  characteristics  of  multistaged,  free-flight,  ballistic  rockets  such  as 
the  DRES  developed  aerial  targets  CRV7/BATS  and  ROBOT-9.  These  equations  of 
motion  form  the  core  of  a  FORTRAN  simulation  software  package  called  BALSIM. 
This  package  allows  for  modeling  of  vehicle  thrust  and  structural  asymmetries,  time- 
varying  mass  and  inertia  characteristics,  variable  wind  conditions,  nonstandard 
atmospheric  conditions,  stage  failures,  and  different  rocket  motor  types.  The  BALSIM 
package  has  been  written  in  IBM  FORTRAN  IV  and  has  been  tested  on  the  IBM  3033 
computer  with  the  H-extended  compiler.  It  is  currently  being  adapted  for  use  with  the 
VAX  11/780  and  Honeywell  DPS-8/70C  computers. 
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1.  INTRODUCTION 

In  1979  DRES  began  development  work  towards  modifying  the  U.S.  army  BATS 
target  for  use  with  the  higher  specific  impulse  CRV7  rocket  motors  under  the  auspices  of 
a  joint  US-Canadian  TTCP  agreement.  This  original  work  has  lead  to  a  number  of 
DRES  initiated  activities  including  the  modification  of  the  target  to  permit  rnultistaging 
and  the  development  of  an  all  CRV7,  multistaged  vehicle  referred  to  as  ROBOT-9 
(Figure  1).  As  well,  support  equipment  has  been  developed  for  target  operation  in 
moderately  heavy  seas,  i.e.  ir  conditions  typical  of  Canadian  coastal  waters  (the  ROBOT 
System  development  history  is  summarized  in  more  detail  in  Reference  1). 

To  support  the  development  of  these  free-flight  targets,  computer  simulation 
programs  were  required  that  predicted  the  dynamic  characteristics  of  the  vehicles.  Of 
particular  importance  were  accurate  predictions  of  basic  performance  parameters  (e.g. 
range  and  flight-time),  wind  effects,  effect  of  nonstandard  atmospheric  conditions  and 
the  dynamic  effects  of  launching  from  a  moving  ship  on  a  finite,  nonzero  length 
launcher. 
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No  existing  software  was  available  to  DRES  that  performed  all  of  these  tasks 
while  conveniently  permitting  some  configuration  variation.  As  a  result,  in  the  period 
September,  1980  to  December,  1981,  a  six  degree- of- freedom  simulation  package  was 
written,  debugged,  and  tested  at  DRES.  This  package  was  applied  to  the  evaluation  of 
the  CRV7/BATS  and  ROBOT-9  performance  characteristics  and  safety-envelopes.  It 
was  coded  in  IBM  FORTRAN  and  has  been  used  on  the  IBM  3033  computer  with  the 
H-extended  compiler.  The  package  is  currently  being  installed  on  a  VAX  11/780 
computer  for  use  with  a  FORTRAN  77  compiler,  and  will  be  adapted  for  use  with  the 
Honeywell  DPS-8/70C  computer.  The  package  has  been  designated  BALSIM. 

It  is  the  intent  of  this  report  to  provide  documentation  of  BALSIM  in  sufficient 
detail  to  permit  users  familiar  with  FORTRAN  to  run  the  program.  Chapter  2  develops 
the  dynamic  model  and  summarizes  its  limitations.  Chapter  3  describes  the  BALSIM 
package  in  general  terms.  Finally,  Volume  2  is  intended  as  an  essentially  self- 
contained  userbook  for  the  package,  and  includes  a  listing  of  all  program  modules. 

2.  DYNAMIC  MODEL 

This  section  summarizes  the  key  features  of  the  dynamic  model  programmed  into 
the  BALSIM  package.  The  basic  six  degree- of- freedom  equations  are  derived  in  the 
following  sections. 

2.1  Fundamental  Assumptions 

Several  overall  simplifying  assumptions  have  been  made  in  the  derivation  of  the 
equations  of  motion.  They  are  valid  for  ballistic  rocket  vehicles  that  have  rigid  structures 
and  relatively  short  ranges,  i.e.  less  than  100  km  (50  nm). 

The  assumptions  are  as  follows: 

1.  The  Earth  is  flat  and  any  Earth -fixed  reference  frame  is  inertial. 

2.  The  vehicle  is  a  rigid  body. 

3.  There  are  no  control  surfaces. 

Assumption  3  may  be  readi'y  relaxed  by  adding  the  appropriate  control  terms  into 
the  equations  of  motion. 
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2.2  Reference  Frames,  Rotation  Matrices  and  Angular  Velocities 

In  the  generai  case  both  an  Earth-fixed  (say  F£)  and  an  inertiai  reference  frame 
(say  F,)  must  be  defined.  Because  of  the  first  simplifying  assumption  of  the  previous 
section,  these  reference  frames  become  identical.  Only  Fr  will  be  used  here.  Thus  let  F, 
be  an  Earth-fixed  inertial  reference  frame  whose  origin  is  at  the  launch  site,  whose  x-axis 
points  along  the  projection  of  the  nominal  launch  trajectory  onto  the  Earth’s  surface  and 
whose  z-axis  is  nominally  downwards  (see  Fig.  1).  Thu  y-axis  follows  from  the  right 
hand  rule. 

A  second,  inertial  Earth- fixed  reference  frame  that  is  useful  >s  the  laureher 
reference  frame  F£.  The  origin  of  F£  is  located  at  the  launch  site  with  the  x-axis  pointing 
in  the  launch  direction  and  the  z-axis  being  nominally  downwards  (see  Fig.  1).  The 
y-axis  follows  from  the  right  hand  rule. 

A  third  Earth- fixed  reference  frame  which  is  occasionally  requited  is  that  of  a 
reference  frame  Fr  placed  at  some  distance  from  the  launch  site.  This  reference  frame 
may  be  used  to  compute  the  rocket  aspect  angle  presentations  (e.g.  from  the  training 
ship).  Since  the  location  and  orientation  of  this  reference  will  depend  on  the  particular 
application,  it  is  defined  only  generally  in  Figure  1. 

Since  the  aerodynamic  forces  are  most  conveniently  expressed  with  respect  to  the 
vehicle,  a  body- fixed  reference  frame  F»  will  also  be  used.  The  origin  of  FB  is  located  at 
the  vehicle  centre- of- mass.  In  vehicles  that  are  axisymmetric,  the  x-axis  is  on  the  ax's  of 
symmetry  and  points  forward  through  the  nose.  Otherwise  the  x-axis  points  in  the 
nominal  launch  direction.  The  z-axis  is  nominally  downward,  while  the  y-axis  follows 
from  the  right  hand  rule.  This  reference  frame  and  some  associated  aerodynamic  angles 
are  shown  in  Figure  2. 

For  cases  where  the  rocket  vehicle  mass  characteristics  are  axisymmetric,  the 
aerodynamic  forces  are  independent  of  the  vehicle’s  roll  attitude,  and  the  thrust  forces 
are  axisymmetric,  the  body- fixed  reference  frame  need  not  spin  with  the  vehicle.  Thus  a 
reference  frame  Ffl'  is  defined  which  is  identical  to  FB  except  that  it  does  not  rotate  with 
the  vehicle  about  the  axis  of  symmetry.  Initially  Ffl  and  Fs'  will  coincide. 

A  third  body-fixed  reference  frame  that  is  useful  in  specifying  the  vehicle’s  mass, 
inertia  and  configuration  characteristics  is  a  reference  frame  FR  whose  origin  is  located 
on  a  nose  datum  plane  on  the  vehicle.  If  the  vehicle  is  axisymmetric,  then  the  origin  of  F« 
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is  on  the  axis  of  symmetry,  and  its  x-axis  points  towards  the  rear  of  the  vehicle  on  the  axis 
of  symmetry.  Otherwise,  the  origin  of  FB  may  be  any  convenient  location  on  the  nose 
datum  plane,  and  the  x-axis  nominally  points  towards  the  rear  of  the  vehicle.  The  z-axis 
is  nominal'y  upward  and  the  y-axis  follows  from  the  right  hand  rule  (see  Figure  2). 

In  the  following,  use  is  made  of  a  number  of  notation  conventions.  In  particular 
if  bBA  denotes  a  rotation  matrix  relating  the  components  of  a  vector  V  expressed  in 
F/i(v^)  to  the  components  of  the  same  vector  in  FB(vs),  then 

v"  =  Lba  Y_a  (2.2,1) 


The  following  definitions,  geometric  relationships,  and  matrices,  will  be 
employed  in  the  presentation  of  the  equations  of  motion. 

1.  The  rotation  matrix  relating  F/  and  FB: 


L  Bl 


COS  6b  COSXPb 

sin<f>aSin0Bcost,'js 

-cos<pBsinyjB 

cos4«sin0Bcos^B 
-t-  sin<j>Bsin*4<B 

=  [I*/,-,] 


cos0Bsinif>B 

sin<}>Bsin0Bsimj;B 
+  COS<}>B  COSIp  B 

cos<{>Bsin0flsinipB 
-  sin<{>BcosifB 


-sin0B  1 

I 

sm<}>Bcos0B  | 


cos<S>Bcos0B 


(2.2,2a) 


(2.2,2b) 


where  <}>„,  0B,  and  ipB  are  the  Euler  angles  defined  by  Etkin  (Reference  3). 

The  rotation  matrix  relating  F;  and  F„'  follows  from  (2.2,2a)  by  substituting  yB  ’ , 
6b'  and  '  for  ipB,  dB  and  §B  respectively,  and  will  be  denoted  LB',. 

2.  The  rotation  matrix  relating  FL  and  F,: 


L  = 


=  lhr„) 


cos  Qgo 

0 

-sin0Bo 

0 

1 

0 

(2.2,3a) 

sin0Bo 

0 

COS0B„ 

J 

(2.2,3b) 

UNCLASSIFIED 


UNCLASSIFIED 


/5 


3.  The  rotation  matrix  relating  FT  and  F,: 


L 


COS  sinter  0 

-  sinter  cosier  G 

0  0  1 


-  Ur/i(] 


(2.2,4a) 


(2.2,4b) 


4.  The  angular  velocity  of  the  vehicle  with  respect  to  F/  written  as  components  in  FB  and 
F 


=  (pfl,  qB,  r b)t 


(2.2,5a) 


Cl) 


B  t  _ 


(0,  qB' ,  r*')1 


(2.2,5b) 


5.  The  angular  rate  cross-product  matrices  (Reference  3)  in  F*  and  F„' : 


0 

-  r* 

q* 

f\j 

y"  = 

r. 

0 

- Pb 

-Qc 

Pa 

o  i 

— 

— — 1 

0 

-rri' 

q*' 

rViK 

o>8  - 

rj' 

0 

0 

a' 

0 

0 

(2.2,6a) 


(2.2,6b) 


6.  The  airspeed  vector  of  the  vehicle  written  as  components  in  and  Fs' : 


V3  =  (Ub,  va,  wfl)r  (2.2,7a) 

V8'  =  (u„  V,  vfB')T  (2.2,7b) 

7.  The  gioundspeed  vector  of  the  vehicle  with  respect  to  F/  written  as  components  in 
Fb,  Ffl'  and  F,:. 


UNCLASSIFIED 


UNCLASSIFIED 


Ye  -  ('Js£,  VBe,  W beY 


(2.2,8a) 


(Ub£,  V%\  W be')t 


(2.2,8b) 


=  (X/,  y,,  z,)7 


(2.2,8c) 


8.  The  aerodynamic  angles  (see  Figure  2): 


=  arctan  (wB/uB) 


(2.2,9) 


where 


=  arctan  (vB/V„) 


—  ^u.  p2  +  v»s2 
=  'ATfl2  +  VB2  +  WB2 


(2.2,10) 


(2.2,11) 

(2.2,12) 


3-* 


Here  a  is  the  angle  of  attack  of  the  x-axis  of  FB,  p  is  the  sideslip  angle  of  the  x-axis 
of  FB,  V  is  the  airspeed,  and  V„  is  the  magnitude  of  the  airspeed  vector  component  along 
the  x-z  plane  of  FB. 

Analogous  angles  to  a  and  p  may  be  written  in  terms  of  FB'  components  by  direct 
substitution  of  Ffl'  quantities  for  FB  quantities. 


9.  The  geometric  relationships  (see  Figure  2): 


uB  =  V  cos  p  cos  a 


(2.2,13a) 


vB  =  Vsin  p 


(2.2,13b) 


wB  =  V  cos/) sin  a 


(2.2,13c) 


10.  The  wind  velocity  with  respect  to  F,  written  as  components  in  F;,  FB,  and  FB' 


W'  =  (W„  W2,  W3)5 


W  =  (uBfi,  vfl<,  wBj;)r 
Wfi'  =  (uv  vB/,  wBg')3 


(2.2,14a) 

(2.2,14b) 

(2.2,14c) 
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II.  The  acceleration  due  to  gravity  written  as  components  in  F,: 


g'  =  (0,0,  %y 

(2.2,15) 

12.  The  aerodynamic  forces  (not  including  thrust  forces)  written  as  components  in  FB 
and  Fb'  : 

A*  =  (XAb<YAb,ZAb)t 

(2.2,16a) 

A*'  =  (XAb,YAb',YAb'V 

(2.2.16b) 

13.  The  aerodynamic  moments  (not  including  thrust  moments)  written  as  components 
in  Fb  and  Fb' : 

MB  =  (LAb,  M^s,  NXfl)T  (2.2,17a) 

Mba'  =  (LAb,MAb',MAb')t  (2.2,17b) 


14.  The  inertia  matrix  of  the  vehicle  with  respect  to  its  centre-of-mass  expressed  in  FB 


(see  Etkin,  Reference  3)  and  Fb'  : 

\B 

k  XX 

-lB 

l*y 

1 

r*  j 

lB  = 

- 1 % 

lB 

kyy 

-\B 

kyi 

(2.2,18a) 

-If, 

-  \B 

Is 

kzz 

\B 

lxx 

0 

0 

lB'  = 

0 

Is' 

kyy 

0 

(2.2,18b) 

0 

0 

Is' 

hy 

15.  The  total  thrust  forces  written  as  components  in  Fb  and  FB' : 


T*  =  (XrB,  Y rB,ZrBy  (2.2, 19a) 

lB'  =  (XrB,  Ytb',  Ytb'Y  (2.2,19b) 
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16.  The  total  tin  ust  moments  is  components  in  FB  and  FB' : 

M?  =  (CrB*  Mt4,  N;b)7  t2-2,20a) 

M ?*  =  (Lra,  MTb\  Mr;)r  (2.2,20b) 

2.3  Nevfton-Euler  Development  of  the  Genera*  Equations  of  Motion 

Newton-Euler  techniques  begin  vvitb  the  fundamental  equations  (Reference  3; 


F 

=  ma 

(2.3,1) 

h 
— > 

=  M 
— > 

(2.3,2) 

a ^  is  the  acceleration  vector  of  the  body  centre -of-mass  reiatr  e  to  an  inertial  reference 
frame,  h  is  the  angular  momentum  of  the  body  about  its  centre- of- mass,  F  is  the 
external  force  vector  acting  a*  the  centre- of- mass  and  M  is  the  external  moment  vector 
about  the  centre- of- mass.  _F  may  be  written 

F  =  mV£  (2.3,3) 

where  VE  is  the  velocity  vector  of  the  vehicle  with  respect  to  F,  and  m  is  mass  of  the 
vehicle.  An  expression  for  h^  follows  from  the  fundamental  relationship 

h  =  /  [  r  x  r ] dm  (2.3.4a) 

^  mass  ^ 

or 

h  =  /  [  r  x  r  +  r  x  ( o)B  x  r )]  dm  (2.3,4b) 

where  is  the  position  vector  of  an  element  of  mass  dm  of  the  body  with  respect  to  its 
centre-of-mass  (see  Figure  3),  o>B  is  the  angular  velocity  vector  of  FB  with  respect  to  F„ 
V  when  applied  tc  a  vector  represents  rate  of  change  with  respect  to  F;  and  ‘o’  when 
applied  to  a  vector  repiesents  rate  of  change  with  respect  to  F*  (see  Reference  4  for  a 
more  thorough  discuss  on  of  vector  differentiation).  Equation  (2  3,4b)  may  be  written 
in  matrix  notation  as  (replacing  coB  x  r  by  -  r  x  a>B  and  dropping  the  subscript  ‘B’ 
cn  co B  for  the  sake  of  brevity)* 

*  Superscripts  cn  matrix  quantities  refer  to  the  reference  frame  ii.  which  *he  components  of  the 
matrix  are  expressed.  Overscore  refers  to  the  matrix  equivalent  of  the  vector  cross-product. 
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But 


hi  =  J  !?*  r*  -  r'«r’>-]dm 

mass 

rB  =  0 


(2.3.5) 

(2.3.6) 


for  a  rigid  body.  Since  cofl  is  a  constant  with  respect  to  the  integration  in  (2.3,5),  it 
follows  that 


h*  =  |ao)fl  (2.3,7) 

where 

Ifl  =  -/  1*1*  dm  (2.3,8) 

mass 

1*  is,  by  convention,  given  by  (2.2,18a). 

The  externally  applied  force  _F  is  made  up  of  an  aerodynamic  component  A ,  a 
thrust  component  T  and  a  gravitational  component  mg  such  that 

F  =  A  +  T  +  mg  (2.3,9) 

Substituting  (2.3,9)  into  (2.3,1),  the  vector  force  equation  becomes 


m^>=  A  +  T  +  mg  (2.3,10) 

The  externally  applied  moment  M  is  made  up  of  an  aerodynamic  component  M* 
and  a  thrust  component  such  that 

M  +  Mr  (2.3,11) 

Substituting  (2.3,11)  into  (2.3,2)  yields 


h.  =  +  Mr  (2.3,12) 

Other  th'*n  the  gravitational  force,  the  dominant  forces  and  moments  acting  on 
the  aircraft  art  due  to  aerodynamic  causes  and  aie  largely  determined  by  its  orientation 
and  configi  ration..  It  is  accordingly  advantageous  to  write  the  matrix  equations  of 
motion  with  respect  to  a  body- fixed  reference  frame.  This  reference  frame  is  chosen  to 
be  FB.  Furthermore,  this  choice  does  not  introduce  any  gravitational  moments  since  the 
origin  of  and  the  centre- of- mass  of  th.  vehicle  coincide. 


UNCLASSIFIED 
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Thus  the  matrix  force  and  moment  equations  become 

m(  VJ  +  w-Vf)  =  AB  +  TB  + 
and 

jts  +  h*  =  M3  +  m 

Substituting  for  hH  from  (2.3,7),  the  moment  equation  becomes 
XBcofl  +  i*or>  +  uBlBus  =  M3  +  Mr 

Writing  out  the  equation  (2.3,13)  in  scalar  form  yields 
m(uB£  +  qfl'vfi£  -  tbvBe)  -  XAt)  +  XTb  +  mglWu 
m(va£  +  rBuflf  -  p*w*£)  =  YAb  +  YTb  +  mglw„ 
m(wB£  +  pB  >b  -  q*ufl£)  ZAd  +  ZTb  +  mglW33 
for  the  force  equations,  end 

I«P*  ~  ~  I«rB  +  i((p,  -  \fyqB  -  \bxzxb  +  I„(r£  -  q«) 

+  (If,  -  &)i»q«  +  Ify rB pe  -  I*q£p,  =  LAj  +  LTjp 

-  I*  p£  +  l&q*  -  IBrB  -  ifypB  +  i*qB  ~  i*  r*  +  If*  (pi  -  r}) 

+  (1«  -  Ifz)rBrB  +  If, Pb Q n  -  1%, q B  =  M»£  +  Mr/i 

-  If.Ps  -  IfzQa  +  IfxfB  -  if,  pfl  -  If,qs  +  I’r*  +  Ify(ql  -  pi) 

+  (lfy  -  I")psqfl  +  If*qcrfl  -  I"  pfi rB  =  N„e  +  Nrs 


(2.3.13) 

(2.3.14) 


(2.3,15) 


XAn  +  Xtb  +  mgl„u 

(2.3,16a) 

YAb  +  Y  tb  +  mglB/?3 

(2.3,16b) 

ZAg  +  Z  TB  +  mglfl/33 

(2.3,16c) 

(3.2,17a) 

(3.2,17b) 

(3.2,17c) 


for  the  moment  equations. 

Kinematic  equations  are  also  required  for  the  linear  and  rotational  position  of  the 
aircraft.  The  linear  position  equations  follow  from 


Vf  =  Lm V| 


(2.3,18) 
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The  rotational  position  equations  are  the  Euler  angle  rate  equations  and  are 
derived  in  Etkin  (Reference  3).  The  resulting  scaiar  kinematic  equations  of  motion  are 
thus  seen  to  be 


x,  =  lB/nuB£  +  1b/2! vB£  4-  *b/3,wB£  (3.2,19a) 

Yj  =  1b/uUb1.-  +  1b/22v3E  +  1b/32wB£  (3.2, 19b) 

Z/  =  1/j/I3Us£  +  1b/jjVB£.  4-  lB,33wB£  (3.2,19c) 

for  linear  position,  and 

<j>B  =  p*  +  qflsin<JiBtan0£  4-  rBcos<|>Btan0B  (2.3,20a) 

0B  =  qBcos<}>B  -  rBsin<}>B  (2.3,20b) 

Wb  =  [qBsin<j>B  4-  rBcos4B]sec8B  (2.3,20c) 

for  angular  position. 


It  should  be  stressed  that  the  variables  (uB£,  vB£,  wBf)  in  equations  (2.3,16), 
(2.3,17)  and  (2.3,19)  are  the  body-axes  components  of  the  vehicle’s  ground  velocity 
vector.  This  is  not  the  same  as  the  equations  developed  in  Reference  2  where  airspeed 
vector  components  in  body-axes  are  used.  Also,  no  assumptions  have  been  made,  up  to 
this  point,  regarding  vehicle  planes  of  symmetry  and  the  symmetry  of  the  thrust  and 
aerodynamic  forces  and  moments.  Finally,  it  should  be  noted  that  no  assumptions  have 
been  made  about  the  mass  and  inertia  characteristics  of  the  vehicle,  i.e.  in  general 

m  ¥=  0 
and 

if,  *  0 

This,  too,  is  an  added  feature  not  present  in  the  equations  developed  in 
Reference  2. 


UNCLASSIFIED 


UNCLASSIFIED 


/ 12 


A  number  of  simplifications  may  be  added  to  these  equations  if  certain  symmetry 
conditions  are  satisfied.  If  the  vehicle  has  mass  symmetry  about  the  xy  and  xz  planes, 
then 

If,  =  If,  =  If,  =  0  (2.3.21) 

If  the  vehicle  mass  characteristics  are  also  axisymmetric,  then  in  addition  to 
(2.3,21)  we  also  have 

If,  =  If,  (2.3,22) 

Applying  assumption  (2.3,21)  to  the  moment  equations  (2.3,17a)  through  to 
(2.3,17c)  results  in  the  simplified  set  of  equations 

pB  =  [-if,pB  -  (If,  -  If,)rBqfl  +  L Ab  +  LrJ/If,  (2.3,23a) 

qB  =  [  -  If,qB  -  (If,  -  If,)p*rB  +  MAb  +  Mr*] /If,  (2.3,23b) 

vB  =  [-if,rB  -(If,  -  Ifjpaq*  +  N„fl  +  NrJ/If,  (2.3,23c) 

Applying  the  axisymmetry  assumption  (2.3,22)  to  those  equations  simplifies 
(2.3,23a)  even  further  by  eliminating  the  (If,  -  If,)rBqB  term,  i.e. 

pB  =  [-if,pB  +  LAb  +  LTa]  /  If,  (2.3,24) 

If  we  now  make  the  assumptions  that  the  thrust  forces  are  axisymmetric,  that  the 
vehicle  mass  characteristics  are  axisymmetric,  and  that  the  vehicle  aerodynamic 
characteristics  are  independent  of  the  vehicle’s  roll  orientation,  then  we  may  take 
advantage  of  the  simplifications  that  will  result  to  the  equations  of  motion  by  expressing 
.hem  in  the  reference  frame  Ffl'  rather  than  in  FB.  Recall  that  in  Section  2.2  we  defined 
F«'  as  being  identical  to  Fa  except  that  it  does  not  rotate  with  the  vehicle  about  the  axis 
of  symmetry.  As  a  result,  if 

b  =  Pb  i  b  +  qB  j  b  +  rB  kB  (2.3.25) 

^  — >  “► 

is  the  angular  velocity  vector  of  Fb”  relative  to  the  inertial  reference  frame  F;  expressed 
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as  components  in  FB\  then  the  latter  condition  implies  that 

p  s'  =  0  (2.3,26) 

for  all  t. 

Also  in  general  we  will  have 

qi  +  qB  (2.3,27a) 

r S,  *  rB  (2.3,27b) 


For  the  purpose  of  generating  the  aerodynamic  forces,  however,  qB  and  rB  may 
be  treated  as  being  interchangeable  with  q„  and  rB  respectively  because  of  the 
assumption  that  the  aerodynamic  forces  are  independent  of  the  roll  orientation. 
However,  the  resulting  aerodynamic  forces  are  now  expressed  in  Fi  axes  rather  than  in 
Fb  axes. 


There  are  also  aerodynamic  forces  that  are  generated  due  to  the  vehicle’s  rolling 
angular  velocity  pB.  The  latter  is  not  identical  to  pi,  and  thus  an  equation  of  motion  for 
pB  will  still  have  to  be  retained. 


Formally  the  equations  of  motion  in  Fi  may  be  obtained  from  the  equations  of 
motion  in  FB  by  setting  pB  =  0  in  all  equations  except  the  pB  equation,  replacing 
(uB£,  vfl£,  w fi£)  with  (ufl£,  vBe,  w Be),  (qB,  rs)  with  (qi,  ri),  (<j>B,  6B,  yjB )  with 
(4>i,  W,  Wb),  (M U„  NAb)  with  (Mifl,  ML),  (Mrs,  Nr,)  with  (Mi„  Mi,),  (X*„  Y„„  ZAb) 
with  ( XAg ,  YAb,  YAb),  (Xr,,  Y Tb,  ZTg)  with  (Xr,,  Yi„  Yi,),  and  applying  the  assumptions 
discussed  previously.  Finally,  the  pB  equation  from  the  FB  equations  is  retained. 


The  resulting  equations  of  motion  are  as  follows: 


i 


m(ufi£  +  qiw^£  -  rB vB£)  =  XAg  +  XTs  +  mgU,,,  (2.3,28a) 

m(vi£  +  rB  ufl£)  =  YkB  +  Yi,  +  mgU/23  (2.3,28b) 

m(wi£  -  qBuS£)  =  Y;,  +  Ytb  +  mgU,33  (2.3,28c) 
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P  B 

= 

[-I2.P.  +  LAg  +  LrJ/I* 

(2.3,29a) 

q  b 

= 

[- i^qi  +  +  My/I' 

(2.3,29b) 

r  B 

= 

l-i&r;  +  m;b  +  My /i* 

(2.3, 29c) 

X/ 

= 

ls/j  [  Ub£  +  lfl/21  Vb£  +  1b/31  W^£ 

(3.2,30a) 

y / 

= 

1w12Ub£  +  lfi/22  VgF  +  1b/32Wb£ 

(3.2,30b) 

Z/ 

= 

1b/13Ub£  +  lfl/23vfi£  +  1b/33Wb£ 

(3.2,30c) 

= 

q^sin^tanfli  +  racos<}>itan0i 

(2.3,31a) 

Bs 

= 

qicos<|>B  -  rising 

(2.3,31b) 

Wb 

= 

[q^sin^  +  rBCOS<f>i]sec0B 

(2.3,31c) 

2.4  Aerodynamic  Model 

The  equations  of  motion  developed  in  the  previous  section  contain  terms 
(e.g.  XB)  that  represent  the  aerodynamic  forces  acting  on  the  vehicle.  In  this  section 
these  terms  are  defined  as  functions  of  the  vehicle’s  state. 

Although  more  sophisticated  techniques  are  available  (see  the  discussion  in 
Reference  2),  for  the  purposes  of  rigid  body  six  degree-of-freedom  simulation,  it  is 
usually  quite  adequate  to  use  a  quasisteady  aerodynamic  model  based  on  Bryan’s 
aerodynamic  derivative  technique. 

No  attempt  has  been  made  here  to  generalize  this  model  so  that  it  applies  equally 
well  to  all  types  of  air  vehicles.  Rather,  its  form  has  been  simplified  so  that  it  is  suitable 
for  use  only  with  free-flight,  ballistic,  rocket-boosted  vehicles. 

The  resulting  model  expressed  as  aerodynamic  force  and  moment  components  in 
Fb  is  summarized  below.  No  attempt  is  made  to  rationalize  this  model  other  than  to  state 
that  its  use  has  resulted  in  predicted  trajectories  that  are  in  good  agreement  with 
measured  flight  characteristics  (see,  e.g.,  Reference  1)  of  CRV7/BATS  and  ROBOT-9 
vehicles. 
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The  aerodynamic  forces  are  specified  by 


XA„  =  -  C0qnS  (2.4,1a) 

T Ag  =  ^ypPQoS  +  QoS  (CyQ/,„  +  Cy0/in)  pseudo  (2.4,1b) 

^Ag  =  CIQ,  Ot  Qo  S  +  Qo  S(C,tt/i„  +  Clofin)  pseudo  (2.43Ic) 

The  various  quantities  in  these  equations  are  defined  in  the  notation.  It  is 
important  to  note  that  qc  is  the  dynamic  pressure  given  by 

qo  =  ‘/:eV2  (2.4,2) 

where  q  is  the  air  density  and  V  is  the  airspeed  given  by 

V  =  (ui  +  vl  +  w I)*'2  (2.4,3) 


Here  (uB,  vB,  ws)  are  the  airspeed  vector  components  expressed  in  FB,  and  in  the 
presence  of  nonzero  wind  conditions  will  not  be  the  same  as  ( uBjr ,  \Be,  wBf?).  Rather,  they 
will  be  related  to  the  wind  velocity  vector  components  in  ¥B  [(uBg,  vBg,  wB(j)]  through  the 
relationships 


Ug  =  U%  -  uflj  (2.4,4a) 

VB  =  v%  -  vBg  (2.4,4b) 


w B  =  w Be  -  \vBg  (2.4,4c) 

The  latter  may  be  obtained  by  considerng  the  fundamental  vector  relationship 


v£  =  V  +  W  (2.4,5) 

i.e.  the  ground  velocity  vector  V  e  equals  the  airspeed  vector  V  plus  the  wind  velocity 
vector  W. 
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CD  is  the  nondimensional  drag  coefficient,  Cyp  is  the  aerodynamic  derivative 
relating  y- force  due  to  sideslip  angle  p,  CtQ  is  the  aerodynamic  derivative  relating  z- force 
due  to  angle  of  attack  a,  and  S  is  a  reference  area  that  is  usually  the  fuselage  cross- 
sectional  area  for  ballistic  vehicles.  CD,  C yp,  and  Cia  may,  in  general,  be  Mach  number 
and  Reynolds  number  dependent,  although  here  they  are  considered  to  be  only  Mach 
number  dependent,  a  and  p  are  given  with  respect  to  the  x-axis  of  Ffl,  and  from 
geometric  considerations  may  be  shown  to  be  (see  Figure  2) 


where 


a 

=  arctan(wB/us) 

(2.4,6a) 

P 

=  arctan(vB/ V„) 

(2.4,6b) 

v„ 

=  (ui  +  w  ly> 

(2.4,7) 

Differential  equations  may  be  obtained  for  a  and  P  by  differentiating  (2.4,6a)  and 
(2.4,6b)  with  respect  to  time,  with  the  results 

a  =  [wfl/V  -  uB  wfl/  V^cos2**  (2.4,8a) 

and 

P  -■  [VXIv*  -  vs(ufluB  +  wBwB)V„]/ (V*zcosp)  (2.4,8b) 

The  last  terms  on  the  right  hand  sides  of  (2.4,1b)  and  (2.4,1c)  are  pseudo  fin 
terms  and  are  included  to  permit  modeling  of  aerodynamic  asymmetries  (e.g.  due  to 
production  tolerances).  They  do  not  include  any  ot  the  effects  produced  by  the  vehicle’s 
nominal  fin  configuration.  The  latter  have  already  been  included  in  CD,  Cyp  and  CIff. 

The  pseudo  fin  terms  are  defined  as  follows: 


C  = 

^yofin 

—  Cia/,„  6ftn  sin$/in 

(2.4,9a) 

tofirt 

Cla/In  6ftn  COS^/m 

(2.4,9b) 

Cy(l//„  = 

(2.4,10a) 

c  — 

^zoi/in 

&fin  C OS  ^ fin 

(2.4,10b) 
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w/in  =  wBcos<t>/in  +  vB  sin^/.n  (2.4,11) 

otft„  =  arctan(w/in/ufl)  (2.4,12) 

Here  CLafin  is  the  lift  slope  of  the  fin,  6f,„  is  the  cant  angle  of  the  fin,  $,in  is  the 
angular  cylindrical  coordinate  of  the  fin,  e/-n  is  a  body-fin  interference  factor  and  afl„  is 
the  angle  of  attack  of  the  fin.  The  fin  geometry  coordinate  system  is  summarized  in 
Figure  4. 

The  aerodynamic  moment  expressions  are  defined  similarly  to  the  aerodynamic 
force  expressions,  as  follows: 


L>tfl  —  Cip  pBqDSbV  (2  V)  +  Cz^/jnd/.nqjoSb 

_  ^ABnom  ^c*  ~ 

Mab  =  ZABnoJxac  -  xc„)  +  Cmfq»SqDb2/(2V) 

XABzct  +  qDS(Czfl./in otf,n  4-  Ct(t/in)p,'Udo(xaC/in  xce) 

NAb  =  -  YAgnom(xttC  -  X")  +  C„rrBSqDb2/ (2V) 

^Ag  y  eg  q  O  ^  (  CyO/in  "h  C'yafw) pseudo  (XOCyrt  Xcs) 


(2.4,13a) 


(2.4,13b) 


(2.4,13c) 


In  these  expressions  YAj3nom  and  ZAgnom  are  given  by  (2.4,1b)  and  (2.4,1c)  without 
the  pseudo  fin  contributions,  (xcg,  ycg,  zcs)  are  the  coordinates  of  the  vehicle  centre-of- 
mass  in  the  vehicle  structural  reference  fra.  ie  F*  (see  Figure  2),  (xoc,  yac,  zac)  are  the 
coordinates  of  the  vehicle  aerodynamic  centre  in  F«,  (x,c/in,  yac/l„,  zaCfin)  are  the 
coordinates  of  the  aerodynamic  centre  of  the  pseudo  fin  in  FR,  b  is  the  reference  length 
(usually  the  fuselage  diameter  for  ballistic  \  ehicles)  and  the  aerodynamic  derivatives  Cip, 
C/(5 /in,  Cmq,  Ciafin,  Cnr,  Cya/.n  are  defined  in  the  notation  list. 


The  aerodynamic  forces  and  moments  may  also  be  written  as  components  in  FB  by 
making  an  identical  set  cf  substitutions  into  (2.4,1)  and  (2,4,13)  as  used  in  converting  the 
equations  of  motion  written  in  Ffc  to  those  written  in  F0' .  It  is  important  to  note  that 
certain  simplifications  result  because  of  the  underlying  assumptions  used  m  developing 
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the  Ffl'  equations,  i.e.  no  pseudo  fin  terms  may  be  included  and  yh  and  z  'B  axes 
characteristics  are  identical. 

For  the  sake  of  brevity,  the  aerodynamic  force  and  moment  expressions  in  Fs  will 
not  be  given  here. 

2.5  Mass  and  Moments  of  Inertia  Models 

The  equations  of  motion  have  been  written  so  fhai  variations  in  the  vehicle’s  mass 
and  inertia  characteristics  (due  to  rocket  motor  propellant  burn)  are  permitted. 
Component  methods  are  used  to  compute  the  total  vehicle  mass  and  moments  of  inertia. 
The  components  considered  are  the  vehicle  airframe,  the  vehicle  payload,  the  vehicle 
rocket  motors  less  propellant  and  the  rocket  motors’  propellant.  Of  these  components, 
only  -he  propellant  charack'istics  are  considered  to  be  variable  with  time.  Finally,  the 
assumption  has  been  made  that  the  payload,  the  rocket  motors,  and  the  propellant  are 
point  masses. 

Under  these  conditions  the  expressions  for  the  vehicle  mass  and  inertia 
characteristics  are  summarized  below.  These  equations  are  given  for  the  reference 
frame  FB.  Position  coordinates  are  with  respect  to  the  vehicle  structural  reference 
frame  Fs.  The  subscripts  used  to  reference  the  different  components  are  as  follows: 

1)  ‘em’  —  airframe  (empty) 

2)  ‘PL’  —  payload 

3)  ‘Me’  —  rocket  motors  less  propellant 

4)  ‘PR’  —  rocket  motor  propellant 

Nm  is  the  total  number  of  locket  motors.  Other  variables  used  are  defined  precisely  in 
the  notation. 

The  expressions  for  the  mass  and  inertia  characteristics  are  given  by 
(AX|  =  X|  -  xcg,  Ay^  -  y^  -  ycg,  and  so  forth) 

Nm 

m  =  au  +  +  Z  [(m*,),  +  (mPfi),]  (2.5,1) 
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em  +  AZrm)  +  mp£,(Ay  Pi  +  A  Zpi) 

Nm  (2.5,2a) 

+  .  Ij  {(mjvfJ.KAy^)?  +  (AzMe)-]  +  (mp*),^*)?  +  (Azp*)?]} 


^ yyem  "F  ^«m(Axem  +  Azem)  +  ITlp£(AXp£  +  AZp£) 

N„  ^  (2.5,2b) 

+  .  Ij  {(mAfe),[(AxAfe)I2  +  (Az^e);]  +  (mpp)J(Axpp)f  +  (Azp*),2]} 


"I"  r^«m(Axem  +  Aytm)  -f-  mp£(AXpt  4-  Ay^£) 

14  (2.5,2c) 

+  ;  Jj  {(^*0, [(AxMe);  +  (Ay.^)?3  +  (mp*),[(Axp*)2  +  (Ayp*),2]} 


Axem  Ay em  nipx,  Akpl  Ay pl 


N„ 

j  —  j  {(^m«)<( Ax*fe)i (Ay^),  4-  (nipp),- ( Axpp),  ( Avpp), } 

(2.5, 2d) 

^x<fn  Azcm  mPL  AxPi:  Azp^ 

N„ 

j  ?  j  ( Axjvfe),  ( Az^),-  +  (mpp),(AxpR),  (Azpp),} 

(2.5, 2e) 

^Yem  AZem  **“  HI p£  A y p^,  A Zp£ 

N„ 

i  =  j  ^mA0* (^y^e)i (Aza*,),  +  (mpp), ( AypR), ( Azpp), } 

(2.5,2f) 

N„ 

fra,»x™  ■*"  mPtxPi  j ? 2  [(m">*)>(xm*).  +  (mpR),(xpR),]}/m 

(2.5,3a) 

N„ 

{m.my.m  +  irwypi  +  .  Zj  ((mmi),(yme),  +  (mpR)1(ypp),]}/m 

(2.5,3b) 

N„ 

{m.myem  +  mPLyPL  +  .1  [(mme),(zme),  +  (mpR),(zpp),]}/m 

(2.5,3c) 
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Because  it  has  been  assumed  that  the  only  mass  changes  are  due  to  propellant 
burn,  in  these  expressions  the  only  time  variable  quantities  will  be  (mPR)„  (xPR|,  yPR_, 
z pit,)-  If  we  further  assume  that  the  propellant  burns  in  such  a  way  that  the  centre-of- 
mass  of  the  propellant  of  a  given  rocket  motor  does  not  change  significantly  (e.g.  as 
would  be  the  case  in  rocket  motors  that  are  not  end  burners),  then  (xPR|,  yPRi,  zPR)  are 
not  time  variable  and  only  (mPR),  need  be  considered.  The  latter  is  related  to  the  specific 
impulse  of  the  rocket  motor  through  the  relationship  (Reference  5) 

(mRR),  =  T,(t)/(IJPig)  (2.5,4) 

where  g  is  the  acceleration  due  to  gravity,  T,(t)  is  the  thrust  of  the  i-th  rocket  motor  as  a 
function  of  time  t,  and  IiPi  is  the  specific  impulse  of  the  i-th  motor. 

Under  these  assumptions  and  with  (2.5,4),  m,  the  moment  of  inertia  time 
derivatives,  and  (xc<,  ycg,  zcg)  may  be  readily  computed.  For  the  sake  of  brevity,  an 
exhaustive  set  of  equations  will  not  be  given.  Typically  we  have 

*  =  ;  J,  [-T,(t)/(4ig)]  (2.5,5) 

f**  —  —  2rnem(Ayfm  yc*  +  Azrmzcg)  —  2mPi(AyyCg  +  AzPLzcg) 

Nm 

+  ilj  {-  2(mAfJ1[(AyMJ1yc,  +  (AzMe),zcJ  +  (mPR),  (2.5,6a) 

[(^Ypr)i  +  (Azpr),  ]  —  2 (mPR),  [( AyPR),  ycg  +  (AzPR),zcs]} 

—  mem(xc*  Aycm  +  Axrm  ycg)  +  mPi(xC4AyPi  +  yc*AxPi) 

Nm 

-  jIi{-(mM,)1[x«I(AyM,)I  +  ycs(AxMJ,]  (mPR),  (2.5,6b) 

( AxPR), ( AyPR),  —  ( mPR),  [xcs ( AyPR),  +  ycs(AxPR),]} 


Xcs  ~  Jj  T.(t)/(LP,)g](xPR),  -  xcgm}m-‘  (2.5,7) 
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The  mass  and  moment  of  inertia  characteristic  formulations  expressed  in  the 
nonrotating  reference  frame  Fs'  follow  from  the  FB  expressions  by  incorporating  the 
simplifying  assumptions  used  for  writing  the  Fc'  equations  of  motion  (see  Section  2.3). 
In  particular,  we  have 

15-  =  I*  =  I  %  =  0  (2.5,8) 

and 

15-  =  I*  (2.5,9) 

Equations  (2.5,8)  and  (2.5,9)  are  just  the  result  of  the  mass  and  inertia 
axisymmetry  assumption  used  in  developing  the  FB'  equations  of  motion. 

For  the  sake  of  brevity,  the  expressions  for  the  Ffl'  mass  and  inertia  characteristics 
will  not  be  given  explicitly. 

2.6  Thrust  Characteristics 

In  Section  2.3  the  equations  of  motion  were  written  in  the  reference  FB  with  the 
thrust  forces  and  moments  written  generally  as  (Xrfi,  YTg,  ZTg )  and  (LrB,  Mrfl,  Nr,) 
respectively.  In  this  section  these  terms  are  examined  in  more  detail. 

The  force  terms  ( XTg ,  YTg,  ZTg)  depend  on  the  time  domain  thrust  characteristics 
and  the  physical  location  and  orientation  of  the  rocket  motors.  This  data  must  be  known 
a  priori  to  the  simulation  and  is  provided  as  input  data  to  the  computer  program  in  the 
form  of  the  thrust  versus  time  look-up  tables.  The  transformations  used  are  summarized 
in  Appendix  2. 

The  thrust  moments  require  a  somewhat  more  detailed  examination.  They  are 
considered  to  consist  of  two  components: 

1)  A  moment  due  to  the  location  and  orientation  of  the  thrust  vector 
relative  to  FB  (see  Figure  4), 

2)  A  moment  induced  due  to  fixed  vanes  or  nozzle  grooves  onto  which 
the  exhaust  jet  impinges. 

Thus  we  have 

Lr„  =  (Lt-b)c*  +  (LrB)n.  (2.6, 1  a) 
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Mrs  -  (Mr,).,  +  (Mr.)*  (2.6,1b) 

Nr.  -  (Nr„)c.  +  (N/fi)„.  (2.6,1c) 

These  characteristics  are  rocket  motor  specific.  It  is  assumed  that  such  data  is 
available  for  the  rocket  motors  used  in  the  simulation.  It  then  follows  that  once  the 
orientation  and  location  of  the  rocket  motor  thrust  vectors  relative  to  the  vehicle  are 
specified,  enough  information  is  available  to  determine  (L/-a,  M,-s,  Nr<J)  as  given  by 
(2.6,1a)  to  (2.6,1c)  (a  detailed  treatment  is  giver,  in  Appendix  2  of  Volume  2). 

An  assumption  that  has  tacitly  been  made  in  this  description  of  the  thrust  effects 
is  that  Coriolis  forces  and  moments  on  the  vehicle  generated  by  the  rocket  motor  exhaust 
are  negligible.  This  need  not  always  be  the  case,  particularly  for  the  moments,  if  the 
exhaust  mass  flow  rate  hi  and  the  exhaust  velocity  vector  relative  to  the  vehicle  V  R  are 
large.  However,  for  vehicles  in  the  class  of  CRV7/BATS  and  ROBOT-9  using  short 
burn  duration  70  mm  (2.75  inch)  rocket  motors,  these  effects  aie  negligible  and  will  not 
be  considered  further  in  this  report. 

2.7  Vehicle  Kinematic  Restrictions  While  on  Launcher 

The  presence  of  the  launcher  during  the  initial  portion  of  the  flight  places  a 
number  of  kinematic  constraints  on  the  vehule’s  motion.  This  section  considers  these 
constraints  for  a  rail  launcher  such  as  was  used  for  CRV7/BATS  and  ROBOT-9  (see 
Reference  1). 

The  basic  geometrical  quantities  are  defined  in  Figure  5  1  he  equations  of  motion 
while  the  vehicle  is  on  the  rail  are  presented  for  the  following  assumptions: 

1)  The  vehxle  is  mechanically  comtrained  from  tipping  backwards  or 
forwards  by  the  guide  T-boit  until  the  T-bolt  cleais  the  launch  rad 
(i.e.  the  vehicle  is  initially  constrained  to  move  along  the  x-axis  of 
reference  frame  F*).  In  the  case  of  CRV7/BAFS  and  ROBOT-9 
there  is  also  the  launcher  cage  constraining  the  vehicle  for  part  of  its 
travel  on  the  launch  rail  (see  Reference  1). 

2)  The  quantity  sG  represents  the  distance  the  vehicle  must  move  in  the 
A-diiecilon  of  r«  beiOie  il.e  guide  T-boh  clears  .he  iaunch  rail.  The 
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bolt  is  assumed  to  be  back  l'ar  enough  on  the  vehicle  so  that  no 
significant  tip-off  may  occur  after  it  is  clear  of  the  rail  and  prior  to  the 
whole  vehicle  coming  clear. 

3)  The  vehicle  may  not  move  backwards  on  the  launcher  rail  (i.e.  ufl  is 
never  less  than  zero). 

Under  these  assumptions  it  follows  that  if  the  distance  that  the  vehicle  centre- of- 
mass  has  travelled  is)  S  less  than  or  equal  to  sc,  then  the  vehicle  is  physically  constrained 
to  move  only  in  the  launch  rail  direction,  i.e.  for  s  <  sG  we  have 


Vb 

= 

qB  =  r«  =  vB  =  \\B  =  0 

(2.7,1a) 

Ub 

> 

0 

(2.7,1b) 

vfl 

= 

Vb(0) 

(2.7,2a) 

W  B 

= 

wB(0) 

(2.7,2b) 

Pb 

= 

P.(0) 

(2.7,2c) 

q  a 

= 

q«(0) 

(2.1, 2d) 

rB 

= 

rB(0) 

(2.7, 2e) 

The  nonzero  cor.didons  (2.7,2a>  to  (2.7,2e)  allow  for  a  nonstationary  launcher, 
i.e.  as  would  be  the  case  for  a  launch  from  a  ship  in  linear  and  angular  motion. 

The  quantity  s  is  defined  precisely  as 

s  ^  \/(x,  -  x,o)2  4-  (y,  -  y,oV  +  z f  (2.7,3) 

where  (x,,  y,,  z,)  are  the  vehicle  centre- of- mass  coordinates  in  F,  and  (x,o,  y/o,  0)  are  the 
centre-of-mass  coordinates  when  the  vehicle  is  at  rest  on  the  launcher  prior  to  first  stage 
ignition. 

For  s  >  sG,  the  governing  equations  are  the  unconstrained  equations  of  motion 
developed  in  Section  2.3. 
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2.8  Wind  Model 

The  aerodynamic  model  presented  in  Section  2.4  includes  the  wind  velocities  in  FB 
(Fa')  as  (uBs,  vBg,  wBg)  [(u»^,  vBg,  Wfls)],  and  has  tacitly  assumed  that  there  is  no  variation 
of  the  wind  velocity  from  one  point  to  another.  This  is  equivalent  to  assuming  that  the 
wind  induced  aerodynamic  loads  are  determined  by  its  velocity  rate  of  change  acting  at 
the  centre-of-mass  of  the  vehicle,  an  assumption  referred  to  as  the  uniform-gust 
approximation  (References  2  and  3).  This  approximation  is  equivalent  to  assuming  that 
the  wind  velocity  spectral  content  significantly  affecting  the  vehicle  response  is  at 
wavelengths  that  are  greater  than  the  significant  vehicle  dimensions  (Reference  3).  This 
assumption  is  reasonable  when  considering  the  rigid  body  dynamic  response  of  flight 
vehicles,  particularly  for  smaller  vehicles  such  as  ROBOT- 9  and  CRV7/BATS. 

The  wind  velocity  vector  components  relative  to  F,  are  most  conveniently 
expressed  as  components  in  F,,  i.e.  (W,,  W2,  W3).  These  components  may  then  be 
related  to  the  wind  velocity  components  in  FB  with  the  rotation  matrix  LBI  as  given  by 
(2.2,2a),  i.e. 


(uv  vv  w Bgy  =  L„(W„  W2,  W3)r  (2.8,1) 

The  simulation  package  has  provisions  for  inputting  (W„  W2,  W3)  as  functions  of 
altitude.  This  allows  modeling  of  wind  velocity  atmospheric  boundary  layer  effects, 
vehicle  encounters  with  jetstream  regions,  and  so  forth.  As  well,  since  meteorological 
winds  aloft  data  is  usually  given  as  a  function  of  altitude,  simulation  of  measured  wind 
conditions  is  facilitated. 

2.9  Atmospheric  Conditions 

Since  the  ROBOT-9  and  CRV7/BATS  vehicles  have  the  capability  to  achieve 
altitudes  well  above  9000  m  (30,000  ft),  an  atmospheric  model  is  required  that  takes  into 
account  variations  in  density  (p),  temperature  (T^),  pressure  (p,,),  and  the  speed  of 
sound  (a)  as  a  function  of  altitude  above  sea  level  (h45i). 

The  models  used  are  based  on  the  U.S.  standard  atmosphere  (1962),  as  is  common 
practice  in  aeronautical  p.igineering,  and  are  valid  within  the  troposphere,  i.e.  for 
fur/.  ^  11,1  GO  m  (36,000  ft)  (see  Reference  6).  They  are  given  by 
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T  a  =  288.15  -  0.0065  hASL  R  (2.9,1) 

?A  =  101300.(1,,/ 288. 1 5) 5  255  (2.9,2) 

a  =  20.0463  VT7  (2.9,3) 

6  =  0.00348454  p„/T„  (2.9,4) 

where 

R  =  vE/l\\ASL  +  r£)  (2.9,5) 


is  in  degrees  Kelvin,  p*  is  in  Pascals,  a  is  in  meters  per  second,  q  is  in  kilograms  per 
meter  cubed,  is  in  meters,  and  r£  is  the  Earih’s  radius  to  the  sea  level  datum, 

re  =  6.3567658  x  106  m  (2.0855531  X  107  ft). 

From  the  factor  R  it  is  also  convenient  to  compute  the  variation  of  the 
acceleration  due  to  gravity  as  a  function  of  altitude,  i.e. 

g  =  goR2  (2.9,6) 

where  g„  =  9.80667  m/s2  (32.1741  f/s2). 

Provision  has  been  made  in  the  simulation  package  to  vary  the  temperature  and 
pressure  (and  thus  the  density)  from  the  standard  values  by  allowing  altitude  dependent 
per  cent  deviations  from  standard  conditions. 

2.10  Aspect  Angle  Equations 

For  target  and  Right  test  applications,  it  is  frequently  necessary  that  the  vehicle’s 
aspect  azimuth  (%A)  and  elevation  (|£x  angles  be  known  with  respect  to  an  observer  at  F7 
(see  Figure  1).  This  section  presents  equations  for  lA  and  4£  in  terms  of  the  location  and 
orientation  of  Fr  relative  to  that  of  F;. 

It  is  assumed  that  the  x-y  planes  of  F,  and  Fr  are  parallel,  i.e.  that  F7  may  be 
rotated  to  Fr  through  a  rotation  ipT  about  the  z-axis  of  F7. 
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Let  the  vector  position  of  Fr  relative  to  F;  be  R  n,  and  that  of  tne  vehicle  centre- 
of-mass  relative  to  the  origin  of  F;  be  R  (see  Figure  6).  It  follows  that  the  vector 
position  of  the  vehicle  relative  to  Fr  is  given  by 

Rr  =  R  -  R„  (2.10,1) 

or  in  matrix  notation 

Rf  =  L„Rf  -  L 77 R ti  (2.10,2) 

where  Lr/  is  the  rotation  matrix  rotating  vector  components  in  F,  to  components  in  FT, 


and  is  given  by  (2.2,4a). 

Equations  (2.10,2)  may  be  written  in  scalar  form  as 
xT  =  (x,  -  xT,)cost +  (yf  -  yr,)sini^r  (2.10,3a) 

yr  =  -  (X/  -  Xti)  sinifiT  +  (y,  -  y„)cosi^r  (2.10,3b) 

Zt  —  Z/  —  z  ti  (2.10,3c) 

From  the  definition  of  and  |£  in  Figure  1 ,  it  follows  that 

4^  =  arctan  (yr/xr)  (2.10,4a) 

|£  =  arctan  (- zT/xr)  (2.10,4b) 


3.  BALSIM  SOFTWARE  DESCRIPTION  —  GENERAL 

The  dynamic  model  described  in  the  previous  chapter  has  been  implemented  in  the 
BALSIM  simulation  package.  AH  coding  was  carried  out  using  IBM  FORTRAN  for  the 
H -extended  compiler.  The  package  has  been  debugged  and  tested  on  the  IBM  3033 
computer,  and  has  been  used  to  predict  the  dynamic  characteristics  of  the  CRV7/BATS 
and  ROBOT-9  vehicles. 
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The  software  is  currently  being  installed  on  VAX1 1/780  and  Honeywell 
DPS-8/70C  computers. 

The  software  consists  of  a  MAIN  program  plus  nine  subroutines  making  up 
approximately  885  FORTRAN  source  statements.  There  are  no  subroutines  or 
functions,  other  than  these,  that  are  not  available  in  standard  FORTRAN  on-line 
libraries. 

The  software  userbook  is  given  in  Appendix  1  of  Volume  2  with  a  source 
language  listing  of  the  package. 

3.1  Software  Capabilities 

The  dynamic  model  implemented  with  the  BALSIM  package  has  already  been 
discussed  in  detail  in  the  previous  chapter.  Its  limitations  will  not  be  considered  further 
here. 


The  package  was  developed  with  the  objective  of  providing  a  convenient  basis  for 
inputting  the  characteristics  of  multistaged  rocket  vehicles  and  predicting  their  dynamic 
rigid  body  characteristics.  By  appropriately  modifying  the  input  data  set,  it  provides  for 

1)  nominal  and  off-nominal  vehicle  mass,  inertia  and  thrust 
characteristics, 

2)  different  motor  types, 

3)  Mach  number  dependent  aerodynamic  characteristics, 

4)  structural  production  tolerances, 

5)  system  failures  (e.g.  stage  and  fin  failures), 

6)  moving  launchers, 

7)  user  specified  initial  conditions, 

8)  user  specified  payload  characteristics, 

9)  tabular  output  in  either  metric  or  English  units,  and 

10)  multiple  case  runs. 

As  well,  with  minor  software  modification,  response  calculations  may  be  stored 
on  disk  for  subsequent  use  with  other  software  (e.g.  plotting  software). 
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3.2  Software  Limitations 

Lx  its  current  form  the  software  is  not  intended  for  use  in  the  following  types  of 
simulations : 

1)  Ballistic  rocket  vehicles  with  control  surfaces. 

2)  Winged  flight  vehicles. 

3)  N^nrigid  vehicles, 

4)  Simulations  where  the  inertial  flat  Earth  approximations  arc  invalid. 

5)  Vehicles  where  the  staging  process  involves  physically  releasing  rocket 
motor  stages. 

Limitations  (1)  and  (5)  may  be  removed  with  relatively  minor  alterations  to  the 
dynamic  model  of  Chapter  2  with  coi  responding  changes  to  the  software. 

3.3  Numerical  Integration  Algorithm 

The  numerical  integration  algorithm  used  to  solve  the  system  of  ordinary 
differential  equations  describing  the  vehicle’s  dynamics  is  a  fixed  step- size,  fourth  order 
Runge-Kutta  method  (sec  Reference  7).  Provision  has  been  made  for  specification 
of  two  step  sizes,  one  for  use  during  racket  motor  burns,  am’  the  other  for  use  during 
coasting  flight.  The  latter  technique  was  found  to  consume: ably  reduce  CPU  time  in 
certain  simulations. 

3.4  SoJiware  Testing  and  Execution  Times 

The  BALSIM  package  has  beer  used  extensively  to  pi  edict  the  performance  and 
dynamic  characteristics  of  the  CRV7/BAfS  and  FOBOT-9  vehicles.  These  predictions 
have  been  used  to  define  the  nominal,  dispeision  and  safe  y-envelope  characteristics  (see 
References  8  and  9)  of  these  vehicles.. 

Flight  test  data  obtained  early  in  the  development  of  CRV7/BATS  (see 
Reference  1)  was  used  to  update  and  validate  the  aerodynamic  model  that  has  been 
employed.  More  recent  comparisons  with  flight  test  data  have  also  proven  to  be 
satisfactory  (also  see  Reference  1). 
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BALSIM  predictions  have  also  been  evaluated  for  consistency  by  comparing 
results  obtained  using  the  equations  of  motion  written  in  F„  with  those  written  in  F*' , 
with  satisfactory  results. 

The  execution  CPU  time  of  the  package  will  depend  on  the  computer  used,  on  the 
step  sizes  chosen,  and  on  the  duration  of  the  flight  time  simulated.  For  the  IBM  2033 
computer,  the  following  CPU  execution  i>«»s  wca  observed  for  simulations  of  the 
ROBOT  9  vehicle  using  a  0.05  second  integration  step  size,  and  a  1.0  second  tabulated 
output  increment: 

1)  Foi  8  cases  averaging  108  simulated  flight  seconds  per  case,  the 
execution  CPU  time  was  343.4  seconds,  yielding  a  0.4  seconds  CPU 
execution  time  per  simulated  flight  second  ratio. 

2)  The  compilation  CPU  time  with  the  H-extended  compiler  was 
25  seconds.; 

3)  The  linkage  editor  CPU  time  was  1 .6  seconds. 

This  completes  the  general  description  of  the  BALSIM  software  package. 
Detailed  user  related  data  is  given  in  Volume  2. 

4.  SUMMARY 

Six  degree-of- freedom,  rigid  body  equations  of  motion  suitable  for  modeling  the 
dynamic  characterises  of  multistaged,  free-flight,  ballistic  rockets  have  been  rigorously 
developed,  and  have  been  implemented  in  a  FORTRAN  software  package  called 
BALSIM.  This  package  allows  for  modeling  of  vehicle  thrust  and  structural 
asymmetries,  time-varying  mass  and  inertia  characteristics,  variable  wind  conditions, 
nonstandard  atmospheric  conditions,  stage  failures,  and  different  rocket  motor  types. 

The  BALSIM  package  has.  been  successfully  used  to  predict  the  performance  and 
dynamic  characteristics  of  the  CRV7/BATS  and  ROBOT- 9  vehicles  both  with  and 
without  moving  launchers.  It  will  be  adapted  for  use  with  the  VAX  11/780  and 
Honeywell  DPS-8  computers  in  the  near  future. 
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